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The Standard Model includes neutrinos as massless particles, but neu- 
trino oscillations showed that neutrinos are not massless. A simple ex- 
tension of adding gauge singlet fermions to the particle spectrum allows 
normal Yukawa mass terms for neutrinos. The smallness of the neutrino 
masses can be well understood within the seesaw mechanism. We analyse 
two cases of the minimal extension of the standard model when one or 
two right-handed fields are added to the three left-handed fields. A second 
Higgs doublet is included in our model. We calculate the one-loop radia- 
tive corrections to the mass parameters which produce mass terms for the 
neutral leptons. In both cases we numerically analyse light neutrino masses 
as functions of the heavy neutrinos masses. Parameters of the model are 
varied to find light neutrino masses that are compatible with experimental 
data of solar Amji and atmospheric Am^ tm neutrino oscillations for normal 
and inverted hierarchy. 
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1. The model 

We extend the Standard Model (SM) by adding a second Higgs doublet 
and right-handed neutrino fields. The Yukawa Lagrangian of the leptons is 
expressed by 

2 

£y = -J2 + H"R A k) D L + H.c. (1.1) 

fc=l 



in a vector and matrix notation, where <3>fc = iT2$t. In expression (1.1) 
ur, and Dl = (y^ £l) T are the vectors of the right-handed charged 
leptons, of the right-handed neutrino singlets, and of the left-handed lepton 
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doublets, respectively, and &k,k = 1,2 are the two Higgs doublets. The 
Yukawa coupling matrices Tk are x ni, while the are ur x n^. 

In this model, spontaneous symmetry breaking of the SM gauge group 
is achieved by the vacuum expectation values ($fc) va c = (0,Vk/V2) , k = 
1,2. By a unitary rotation of the Higgs doublets we can achieve ('I^vac = 
vl\j2 > and (^vac = with v ~ 246 GeV. The charged-lepton mass 
matrix Mi and the Dirac neutrino mass matrix Mo are 

M e = ^=T 1 and M D = ^=A ll (1.2) 

with the assumption that Mi = diag (m e , m^, m T ). The hat indicates that 
Mi is a diagonal matrix. The mass terms for the neutrinos can be written 
in a compact form with a (til + tir) x (m + ur) symmetric mass matrix 

M,= ( ° 5 V (1-8) 



. M D Mr 

M v can be diagonalized as 

U T M U U = m = diag (mi, m 2 , . . . , (1-4) 

where the rrij are real and non-negative. In order to implement the seesaw 
mechanism [H [2] we assume that the elements of Mb are of order mo 
and those of Mr are of order ijir, with m/j <C tti/j. Then the neutrino 
masses 77tj with i = 1,2, . . . ,til are of order m D /mR, while those with 
i = ni, + 1, . . . , hl + ur are of order uir. It is useful to decompose the 
(til + ttr) x (til + ur) unitary matrix U as U = (Ul,U r ) T , where the 
submatrix Ul is ni x [ni+nR) and the submatrix JTr is ur x (77^+77^) OH]. 
With these submatrices, the left- and right-handed neutrinos are written as 
linear superpositions of the 77 l + ur physical Major ana neutrino fields xf- 
vl = UlPlX an d vr = UrPrx, where Pl and Pr are the projectors of 
chirality. 

It is possible to express the couplings of the model in terms of the mass 
eigenfields, where three neutral particles namely, the Z boson, the neutral 
Goldstone boson G° and the Higgs bosons H® couple to neutrinos. The full 
formalism for the scalar sector of the multi-Higgs-doublet SM is given in 
Refs. |310| . 

Once the one-loop corrections are taken into account the neutral fermion 
mass matrix is given by 

(i) = ( 5M L M T D + 5M T D \^( 5M L M T D \ 

v \ M D + 5M D Mr + 5M r )~\ M D Mr r 1 ' > 
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where the O3X3 matrix appearing at tree level (1.3) is replaced by the contri- 
bution 5 Ml. This correction is a symmetric matrix, it dominates among all 
the sub- matrices of corrections. Neglecting the sub-dominant pieces in ( 1.5 ), 
one-loop corrections to the neutrino masses originate via the self-energy 

function Sf (0) = E S L (Z) (0) +Z S L {G °\o) +Z S L ( H °\o), where the Z S L {Z ' G °' H °\o) 
functions arise from the self-energy Feynman diagrams and are evaluated at 
zero external momentum squared. Each diagram contains a divergent piece 
but when summing up the three contributions the result turns out to be 
finite. 

The final expression for one-loop corrections is given by [5] 



+ ^ M - u 'H^- i TH^h MB ' (i - 6> 

with A;, = 6fcAfc, where b are two-dimensional complex unit vectors and 
the sum runs over all neutral physical Higgses H®. 



2. Case ur = 1 

First we consider the minimal extension of the standard model by adding 
only one right-handed field vr to the three left-handed fields contained in 
vl- 

For this case we use the parametrization Aj = (^/2mL)/v) of, where 
af = (0, 0, l) an d a5T = (0, 1, e^. Diagonalization of the symmetric mass 



matrix M v (1.3) in block form is 



U T M U U = U T ( ° 3x i T "tf 1 )U = ( Ml ? ) . (2.1) 
1 m D a( M R J V M h 1 



The non zero masses in Mi and Mh are determined analytically by fin- 
ding eigenvalues of the hermitian matrix M U M„. These eigenvalues are the 
squares of the masses of the neutrinos Mi = diag(0, 0,m;) and Mh = mh- 
Solutions m 2 D = rrihmi and m\ = (mh — mi) 2 « m\ correspond to the 
seesaw mechanism. 
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Fig. 1. Calculated masses of two light neutrinos as a function of the heavy neutrino 
mass m/j. Solid lines show the boundaries of allowed neutrino mass ranges when the 
model parameters are constrained by the experimental data on neutrino oscillations 
with # a tm — 45°. The allowed values of m; 1 and m; 2 form bands, their scattered 
values are shown separately in the middle plots. 



The diagonalization matrix U for the tree level is constructed from a 
rotation matrix and a diagonal matrix of phases [/tree = C^34(/3)C/j, where 
the angle j3 is determined by the masses mi and m^. 

For the calculation of radiative corrections we use the following set of 
orthogonal complex vectors: bz = (i,0), b\ = (1,0), 62 = (0,z) and 63 = 
(0,1). Diagonalization of the mass matrix after calculation of one-loop 
corrections is performed with a unitary matrix Z7i oop = C/ egv C/ ¥ ,((^i, if2, ^3), 

where U egv is an eigenmatrix of and Uu, is a phase matrix. The 

second light neutrino obtains its mass from radiative corrections. The third 
light neutrino remains massless. 

The masses of the neutrinos are restricted by experimental data of solar 
and atmospheric neutrino oscillations [6j and by cosmological observations. 
The numerical analysis shows that we can reach the allowed neutrino mass 
ranges for a heavy singlet with the mass close to 10 4 GeV and with the 
angle of oscillations fixed to 6> atm = 45°, see Fig. [T] The free parameters 
m H o , m H o , and 4> are restricted by the parametrization used and by os- 
cillation data. Figure [2] illustrates the allowed values of Higgs masses for 
different values of the heavy singlet. 
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Fig. 2. The values of the free parameters m H o and m H o as functions of the heaviest 
right-handed neutrino mass rrih, for the case = 1. The mass of the SM Higgs 
boson is fixed to m H o — 125 GeV and the angle of oscillations is atm = 45°. 



3. Case tir = 2 

When we add two singlet fields vr to the three left-handed fields vl, the 
radiative corrections give masses to all three light neutrinos. 

Now we parametrize Aj = ^ (m^af^D^) with |ai| = 1, = 1, 



\a,2\ = 1 and l^l = 1. Diagonalizing the symmetric mass matrix M u (1.3) 
in block form we write: 

rr, I ° 3x3 m D2 a m Dl b \ / /> n \ 

U MyU — U I raa g „ b \v=(« I ). (3.!) 

The non zero masses in Mi and are determined by the seesaw mecha- 
nism: m 2 D . ~ m^ra^ and m 2 R _ ~ m?., i = 1,2. Here we use mi > > 
ordering of masses. The third light neutrino is massless at tree level. 

The diagonalization matrix for tree level i/tree = U^fU^ifii) is composed 

of an eigenmatrix of M v Ml and a diagonal phase matrix, respectively. 

For calculation of radiative corrections we use the same set of ortho- 
gonal complex vectors bi as in the first case. Diagonalization of the mass 
matrix including the one-loop corrections is performed with a unitary ma- 
trix t/ioop = Ue^Uiplipi), where lllg° p is the eigenmatrix of My and 
Utp is a phase matrix. 
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A broader description of the case n# = 2 and graphical illustrations of 
the obtained light neutrino mass spectra is given in Ref. |7]- Both normal 
and inverted neutrino mass orderings are considered. 

4. Conclusions 

For the case ur = 1 we can match the differences of the calculated light 
neutrino masses to AmQ and Am^ tm with the mass of a heavy singlet close 
to 10 4 GeV. The parametrization used for this case and restrictions from the 
neutrino oscillation data limit the values of free parameters. Only normal 
ordering of neutrino masses is possible. 

In the case tir = 2 we obtain three non vanishing masses of light neutri- 
nos for normal and inverted hierarchies. The numerical analysis [7] shows 
that the values of light neutrino masses (especially of the lightest mass) de- 
pend on the choice of the heavy neutrino masses. The radiative corrections 
generate the lightest neutrino mass and have a big impact on the second 
lightest neutrino mass. 
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